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ABSTRACT 

Stardust will be the first mission to return samples from 
beyond the Earth-Moon system. The sample return cap- 
sule, which is passively controlled during the fastest Earth 
entry ever, will land by parachute in Utah. The present 
study analyzes the entry, descent, and landing of the re- 
turning sample capsule. The effects of two aerodynamic 
instabilities are revealed (one in the high altitude free 
molecular regime and the other in the transonic/subson- 
ic flow regime). These instabilities could lead to unac- 
ceptably large excursions in the angle-of-attack near peak 
heating and main parachute deployment, respectively. To 
reduce the excursions resulting from the high altitude 
instability, the entry spin rate of the capsule is increased. 
To stabilize the excursions from the transonic/subsonic 
instability, a drogue chute with deployment triggered by 
an accelerometer and timer is added prior to main para- 
chute deployment. A Monte Carlo dispersion analysis of 
the modified entry (from which the impact of off-nomi- 
nal conditions during the entry is ascertained) shows that 
the capsule attitude excursions near peak heating and 
drogue chute deployment are within Stardust program 
limits. Additionally, the size of the resulting 3-0 landing 
ellipse is 83.5 km in downrange by 29.2 km in cross- 
range, which is within the Utah Test and Training Range 
boundaries. 

NOMENCLATURE 


CFD 

Computational Fluid Dynamics 

c-g. 

Center-of-Gravity 

DOF 

Degree-of-Freedom 

DSMC 

Direct Simulation Monte Carlo 

LAURA 

Langley Aerothermodynamic Upwind 
Relaxation Algorithm 
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NASA 

National Aeronautics and Space 
Administration 

NNE 

North-Northeast 

PICA 

Phenolic Impregnated Carbon Ablator 

POST 

Program to Optimize Simulated Trajectories 

SRC 

Sample Return Capsule 

ssw 

South-Southwest 

TLNS3D 

Thin Layer Navier Stokes Three- 
Dimensional 

UTTR 

Utah Test and Training Range 

0^ 

Total angle-of-attack (angle between the 
spin axis and velocity vector), deg 

y 

Flight-path angle, deg 


INTRODUCTION 


The fourth of NASA’s Discovery class missions is a 
comet sample return mission known as Stardust. It will 
be the first mission to return samples from beyond the 
Earth-Moon system. The spacecraft is scheduled to be 
launched in February of 1 999 for encounter with the com- 
et Wild-2 in 2004. Stardust will come within 100 km of 
the comet nucleus and deploy a sample tray to collect 
cometary and interstellar dust particles (Fig. la). Upon 
Earth return in January 2006, the entry capsule (Fig. lb), 
containing the comet samples, will be released from the 
spacecraft and land by parachute at the Utah Test and 
Training Range (UTTR). The entry velocity will be the 
highest of any Earth-returning mission (relative velocity 



Fig. la. Stardust spacecraft flight configuration. 




Fig. lb. Stardust sample return capsule configuration. 

about 12.6 km/s). A new heat shield made of PICA (Phe- 
nolic Impregnated Carbon Ablator) will be used to pro- 
tect the Sample Return Capsule (SRC) from the intense 
heat of reentry. 1 

Four hours prior to entry, the SRC will be spun-up 
and separated from the main bus. The SRC has no ac- 
tive control system, so the spin-up is required to main- 
tain its entry attitude (nominal 0° angle-of-attack) dur- 
ing coast. Throughout the atmospheric entry, the pas- 
sive SRC will rely solely on aerodynamic stability for 
performing a controlled descent through all aerodynam- 
ic flight regimes: hypersonic-rarefied, hypersonic-tran- 
sitional, hypersonic-continuum, supersonic, transonic, 
and subsonic. The SRC must possess sufficient aerody- 
namic stability to overcome the gyroscopic (spin) sta- 
bility in order to minimize any angle-of-attack excur- 
sions during the severe heating environment. Addition- 
ally, this stability must persist through the transonic and 
subsonic regimes to maintain a controlled attitude at 
parachute deployment. 

The objective of this study is to analyze the entry, 
descent, and landing of the returning sample capsule. 
This analysis consists of performing a trajectory simu- 
lation of the entire entry (from bus separation to land- 
ing) to predict the descent attitude and landing conditions. 
In addition, a Monte Carlo dispersion analysis is per- 
formed to ascertain the impact of off-nominal conditions 
which may arise during the entry to determine the robust- 
ness of the Stardust SRC design. Specifically, the SRC 
attitude near peak heating and parachute deployment is 
of interest, along with the landing footprint ellipse. 

The SRC is restricted to land within the UTTR site. 
For mission success, a high-fidelity aerodynamic data- 
base is essential in order to accurately predict the land- 
ing location, as well as the attitude of the SRC at critical 
phases (e.g., peak heating and parachute deployment) 
during the entry. In this paper, the aerodynamics utilized 


in the entry simulation is mentioned first, followed by a 
description of the nominal entry sequence of the SRC. 
Finally, the results of the Monte Carlo entry dispersion 
analysis are presented. 

ANALYSIS 


Aerodynamics 

The aerodynamic database utilized for the SRC in 
the flight simulation studies is constructed from a com- 
bination of computational fluid dynamics (CFD) calcu- 
lations, and wind tunnel and historical flight data as de- 
scribed by Mitcheltree, et. al. 2 This large variety of sources 
for the aerodynamics is required because the SRC 
traverses many different flow regimes (hypersonic-rar- 
efied, hypersonic-transitional, hypersonic-continuum, 
supersonic, transonic, and subsonic) during its entry. At 
the outer reaches of the atmosphere, free molecular flow 
aerodynamics are employed. In the rarefied flow regime, 
Direct Simulation Monte Carlo (DSMC) simulations up 
to angles-of-attack of 30 degrees are used to define bridg- 
ing functions for the aerodynamic coefficients. In the 
hypersonic-continuum regime, a matrix of 22 solutions 
from the computational fluid dynamics code LAURA 
(Langley Aerothermodynamic Upwind Relaxation Al- 
gorithm) 3 describe the aerodynamics for Mach numbers 
between 42 and 7.15. At supersonic and transonic speeds, 
the aerodynamics are taken from two sets of existing 
wind tunnel data, augmented by CFD results from the 
TLNS3D (Thin Layer Navier Stokes Three-Dimension- 
al) code. 4 Subsonic aerodynamics are defined by a com- 
bination of static wind tunnel measurements and dynam- 
ic free flight measurements. 5 These sources are blended 
to form a comprehensive database which describes the 
aerodynamics of the SRC for the expected flight condi- 
tions. Figure 2 shows the range of application of the var- 
ious aerodynamic sources mentioned above. The aerody- 
namic characteristics of the SRC are described in detail 
in Ref. 2. 
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Fig. 2 Stardust SRC aerodynamic database. 


2 



Trajectory Simulation 


Table 2. Atmospheric Mission Uncertainties 


The trajectory analysis is performed using the six- 
and three-DOF (degree-of-freedom) versions of the Pro- 
gram to Optimize Simulated Trajectories (POST). 6 This 
program has been utilized previously for similar appli- 
cations. 7,8,9 The three-DOF program (which integrates 
the translation equations of motion) is used from bus 
separation to atmospheric interface. The six-DOF ver- 
sion of POST (which integrates the translational and ro- 
tational equations of motion) is used from atmospheric 
interface until parachute deployment. The three-DOF 
program is used again from parachute deployment to 
landing. The trajectory simulation includes Earth atmo- 
spheric (GRAM-95) 10 and gravitational models, capsule 
separation and non-instantaneous parachute deployment 
models, and capsule aerodynamics and mass properties. 
The validity of the present approach has been demon- 
strated recently through comparisons between the Mars 
Pathfinder pre-flight predictions of the flight dynamics 
and the flight data. 11 

During the entry, off-nominal conditions may arise 
which affect the descent profile. These off-nominal con- 
ditions can originate from numerous sources, such as 
capsule mass property measurement uncertainties, sep- 
aration attitude and attitude rate uncertainties, and lim- 
ited knowledge of the flight-day atmospheric properties 
(density, pressure, and winds). Additionally, computa- 
tional uncertainty with the aerodynamic analysis and un- 
certainties with parachute deployment are contributing 
sources of uncertainty. In this analysis, an attempt is 

Table L Exo- Atmospheric Mission Uncertainties. 


Mass Properties 3-c Variance 

Mass ±0 5 kg 

eg position along spin axis ±0.254 cm 

eg position off spin axis ±0.254 cm 

Major moment of inertia (I**, lyy, l zz ) ±20% 

Cross products of inertia (l xy , l^, l yz ) ±0.015 kg-m 2 


Post-Separation State Vector 

Position 1 correlated with coveriance 

Velocity J matrix corresponding to At) =±0.075 deg 


Pitch attitude ±2.0 deg 

Yaw rate ±6.0 deg/s 

Roll rate +4 rpm, -2 rpm 

Separation 

Spring induced velocity: 

radial velocity ±0.0482 m/s 

cross-track velocity ±0.0482 m/s 

in-track velocity ±0.04 m/s 


3-c 

Aerodynamic Variance 

Free molecular aerodynamics, CA ±10% 

CN, CY ±8% 

Cm, Cw ±12% 

Hypersonic continuum aerodynamics, CA ±4% 

CN, CY ±8% 

Cm, Cw ±10% 

Supersonic continuum aerodynamics, CA ±10% 

CN, CY ±5% 

Cm, Cw ±8% 


Subsonic continuum aerodynamics, CA ±5% 

Hypersonic dynamic stability coefficients, 

Cmq, Cnr... ±0.15 

Supersonic dynamic stability coefficients, 

Cmq, Cnr... ±0.15 


Atmosphere 

Pressure, density, winds: GRAM-95 model 3-o 

scale factor 


Other 

Ablation mass ±1 0% 

* Drogue chute g-switch ±1 0% 

* Drogue chute deployment timer ±1% 

* Drogue chute aerodynamics, CA ±1 0% 

* Main chute deployment timer ±1% 

* Main chute aerodynamics, CA ±1 5% 


'Uncertainty sampled using uniform distribution 

made to conservatively quantify and model the degree 
of uncertainty in each mission parameter. For this mis- 
sion, 41 potential uncertainties were identified. These un- 
certainties are grouped into two categories (exo-atmo- 
spheric and atmospheric) and are listed in Tables 1 and 
2, respectively, along with the corresponding 3-o vari- 
ances. For most of the parameters, a Gaussian distribu- 
tion is sampled. However, for the center-of-gravity (c.g.) 
offset quadrant and parachute deployment parameters (g- 
switch, timers, and aerodynamics), uniform distributions 
are utilized to model their operating performance. 

As will be shown in the results, the successful re- 
turn of the cometary samples by the Stardust SRC de- 
pends heavily on the validity of the Monte Carlo analy- 
sis more than any previous mission. Increased reliance 
on entry simulations for mission success places consid- 
erable importance on selecting appropriate uncertainties. 
As confidence increases in the analysis accuracy, cheaper 
and/or higher performance entry systems can be select- 
ed for future missions. 
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RESULTS AND DISCUSSION 
Nominal Mission 

Original Entry Sequence 

In the original nominal Stardust entry sequence, the 
SRC enters the atmosphere with a spin rate of 5 revolu- 
tions per minute (rpm). The spin rate maintains entry 
attitude (nominal 0° angle-of-attack) until atmospheric 
interface (since the SRC possesses no active control sys- 
tem). As the SRC descends, it must rely solely on atmo- 
spheric stability in all flow regimes to minimize any 
angle-of-attack excursions until main parachute deploy- 
ment at Mach 0. 16. However, the SRC was found to be 
statically unstable in the free molecular flow regime due 
to its aft center-of-gravity location (0.283 m or 0.351 
body diameters back from the nose). 2 The six-DOF anal- 
ysis reveals that the pitch rate induced by the instability 
during the free molecular regime carries into the transi- 
tional region where high angles-of-attack are produced. 
This static instability causes the SRC to pitch-up to a 
total angle-of-attack (Oj) above 70 deg as seen in Fig. 3. 
As the SRC descends into the continuum regime (where 
it is statically stable) the angle-of-attack damps out and 
reduces to small values. 



Attitude, km 


Fig. 3. Original nominal Stardust attitude profile . 

High angles-of-attack early in the entry are a con- 
cern, since they can lead to angles-of-attack greater than 
10 deg at peak heating (which is a Stardust program lim- 
it). Angles-of-attack greater than 10 deg increase after- 
body heating near the shoulder regions and can damage 
the afterbody thermal protection system. Moreover, high 
angles-of-attack are worrisome because the SRC is sta- 
ble flying backwards. Off-nominal attitude and attitude 
rate conditions at atmospheric interface could result in a 
backward entry. In addition to the high altitude static 
instability, a low altitude dynamic instability in the tran- 


sonic and subsonic flow regimes also exists, again caused 
by the aft center-of-gravity location. This transonic/sub- 
sonic dynamic instability could induce a tumbling mo- 
tion prior to main parachute deployment, raising con- 
cerns for a successful deployment. Reference 2 describes 
the SRC’s stability/instability in the various regimes in 
more detail. 

Modified Entry Sequence 

The original nominal entry sequence is unaccept- 
able since off-nominal conditions could result in mis- 
sion failure. Modification of either the entry sequence 
or the SRC is required to improve the probability of mis- 
sion success. 

The high altitude and transonic/subsonic instabili- 
ties could be eliminated by moving the c.g. of the SRC 
forward to 0.26 body diameters back from the nose. 2 
Since the SRC has only three major components (fore- 
body heat shield, sample tray, and parachute canister), 
movement of the c.g. via re-packaging of the SRC is 
difficult. The size and mass of the sample tray preclude 
large movements in the c.g. Ballast could be added to 
the nose of the SRC to move the c.g. forward. However, 
a prohibitive large amount (22.4 kg) of ballast is required 
to move the c.g. sufficiently forward to remove the in- 
stabilities. 

Augmentation of the SRC’s stability is deemed nec- 
essary to eliminate the large angle-of-attack excursions. 
Several concepts were considered. For example, adding 
an aft-skirt would provide a restoring torque in the free- 
molecular regime stabilizing the SRC. 12 However, these 
devices, once serving their purpose, must be discarded 
in order to avoid destabilizing the lower flight regimes. 
To avoid such complications, the solution selected to 
address the high altitude instability is to increase the spin 
rate of the SRC upon entry. The higher spin rate, al- 
though not eliminating the instability, increases the gy- 
roscopic stability of the SRC sufficiently to retard the 
effects of the free molecular static instability. However, 
if the entry spin rate is too large, the gyroscopic stability 
could overwhelm the aerodynamic stability in the con- 
tinuum regime. This would lead to large angles-of-at- 
tack during peak heating. After a detailed investigation 
performing numerous six-DOF entry analyses for a va- 
riety of spin rates, an entry spin rate of 16 rpm is select- 
ed. This spin rate adequately reduces the high altitude 
angle-of-attack excursions, yet avoids any attitude con- 
cerns during peak heating. Additionally, a 16 rpm spin 
rate affords sufficient margin (in the angle-of-attack ex- 
cursions) to accommodate off-nominal conditions which 
may be present during the entry (as confirmed by the 
Monte Carlo analysis presented later). 
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The transonic/subsonic instability is addressed by 
deploying a supersonic drogue chute prior to main para- 
chute deployment. The drogue chute serves as a stabi- 
lizing mechanism for the SRC until main parachute de- 
ployment. The drogue chute size and deployment Mach 
number are constrained by the need to prevent exces- 
sive drift, which could lead to a landing footprint be- 
yond the proposed UTTR site. However, the drogue chute 
size must provide sufficient area to stabilize the SRC. 
Furthermore, the deployment Mach number must be 
outside the dynamic instability region near transonic 
speeds to avoid the possibility of large angles-of-attack. 
From spin tunnel tests, the drogue chute is sized to 0.828 
m in diameter to provide ample area for stabilizing the 
SRC. 5 Numerous six- and three-DOF analyses of the 
entry reveal that drogue deployment at Mach 1.4 avoids 
excessive drift concerns (as confirmed by the Monte 
Carlo results presented later). 

Adoption of these changes into the mission required 
modification of the entire terminal descent procedure of 
the entry. A new deployment algorithm, consisting of a 
g-switch and two timers, is utilized for deployment of 
the drogue and main parachutes. Previously, only a baro- 
switch was needed for deploying the main parachute 
(diameter = 8.2 m). Figure 4 shows the modified nomi- 
nal entry profile, with the terminal descent sequence 
highlighted. The g-switch is triggered after sensing 3 g's 
(decelerating side) for 0.5 seconds, at which point, the 
drogue timer is initiated. After 1 5.04 seconds, the drogue 
chute is deployed, initiating the main timer. After 350.6 
seconds, the main parachute is deployed. This new nom- 
inal entry sequence is sufficiently robust to accommo- 
date off-nominal conditions during the entry (as con- 
firmed by the Monte Carlo analysis presented below). 

Trajectory calculations are repeated for the modi- 
fied entry profile using the most current mass properties 



Fig. 4. Modified nominal Stardust entry sequence. 


Table 3. Nominal Mass Properties of the SRC. 


Mass, kg 46.0 

Center of gravity, m 

Along spin axis (x-direction, from nose) 0.2831 

Off spin axis (y-direction) 0.000396 

Off spin axis (z-direction) -0.00271 5 

l xx , kg-m 2 (spin axis) 2.163 

lyy, kg-m 2 1.595 

l zz , kg-m 2 1.4991 

l xy , kg-m 2 0.00181 

l xz , kg-m 2 0.00221 

l yz , kg-m 2 0.00437 


for the SRC (Table 3). The flight characteristics of the 
modified nominal are shown in Figs. 5 and 6. During 
the entry, the SRC aerodynamically decelerates from 12.6 
km/s to subsonic speeds. The maximum deceleration 
experienced by the SRC during the descent is 32.9 g’s. 

Recall that the SRC is still statically unstable in the 
free molecular regime. The higher spin rate only delays 
the effect of the static instability so that the SRC can 
traverse the transitional regime to the stable continuum 
regime without experiencing a large increase in the total 
angle-of-attack. As seen in Fig. 5, the total angle-of-at- 
tack pitches up to approximately 7 deg in the transition- 
al regime before reducing to less than 2 deg near peak 
heating (which occurs around Mach 35.3). Reference 13 
describes the heating environment encountered during 
the entry. 

As the SRC descends, the static margin decreases 
near Mach 12 to produce a new trim point. Consequent- 
ly, since the SRC has a non-zero c.g. off-set from the 



Fig. 5. Modified nominal Stardust attitude profile. 
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(a) Entry sequence. 


Monte Carlo Dispersion Analysis 

Independent Uncertainty Effects 

Before a combination of off-nominal conditions are 
examined, a sensitivity analysis is first performed to iden- 
tify the mission uncertainties which have the greatest 
impact on the overall landing footprint. Each of the 41 
mission uncertainties are varied independently at their 
respective ±3-0 (maximum/minimum) variance. Figure 
7 shows the resulting downrange obtained from the larg- 
est contributors to the overall landing footprint. Those 
mission uncertainties which are not depicted lead to 
downrange dispersions less than 0.5 km. 



(b) Parachute deployment sequence. 
Fig . 6. Modified nominal mission profile. 


30 r 


Mission Uncertainty 

1 Initial state vector 

2 East-West wind 

3 Density 

4 North- Sooth wind 

5 Radial separation velocity (±0.02 m/s) 

6 Vertical wind 

7 C^: hypersonic-continuum (±4%) 

8 C^: supersonic-continuum (±10%) 

9 Initial mass (±0.5 kg) 



3 4 5 6 7 

Mission uncertainty 


Fig. 7. Significant contributors to the total range 
dispersion. (3-<J variance shown in parenthesis.) 


spin axis, an increase in is observed from a mean (X^ 
of approximately 1 .5 deg near Mach 12 to approximate- 
ly a mean GCj of 2.5 deg near Mach 2. In transitioning to 
a new trim point, attitude rates induce an overshoot in 
(Xp (peaking around Mach 8) before receding around 
Mach 2. As the SRC approaches transonic speeds, the 
dynamic instability drives another increase in until 
drogue chute deployment. 

Beginning at Mach 1.4 (approximately 34 km alti- 
tude), the terminal descent phase of the entry begins, 
which slows the SRC down to approximately 4 m/s pri- 
or to landing. Figure 6b shows the nominal altitudes of 
the drogue and main parachute deployments. 

Note, the mass properties for the SRC continue to 
fluctuate as its design matures. However, the current 
properties listed in Table 3 are representative of the fi- 
nal configuration. 


The mission uncertainties shown in Fig. 7 can be 
grouped into two categories: large contributors (mission 
uncertainties 1-4) and small contributors (mission un- 
certainties 5-9). The first group, containing initial state 
vector and atmospheric wind and density uncertainties, 
contribute on the order of 20-25 km each to the landing 
footprint size. Again, since the atmospheric winds have 
a significant impact on the downrange due to parachute 
drift, the selection of an appropriate drogue chute size 
and deployment Mach number is critical. The second 
group, containing uncertainties in initial mass, bus sep- 
aration velocity, and aerodynamic drag, produce down- 
range dispersions of approximately 1-10 km each. Ta- 
ble 4 summarizes the independent ±3-0 dispersion re- 
sults. Note that uncertainties in the aerodynamics asso- 
ciated with the mission have a minimal impact on the 
overall landing footprint. 

Multiple Uncertainty Effects 

To determine the robustness of the Stardust SRC 
entry profile, off-nominal conditions are simulated to 
address uncertainties which may arise during the descent. 
The impact of multiple uncertainties occurring simulta- 
neously is ascertained by performing a Monte Carlo dis- 
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Table 4. Major Contributors to Total Downrange 
Dispersion * 



Dispersion 

Dispersion 


with +3-o 

with -3-o 


uncertainty, 

uncertainty, 


km 

km 

1 . State vector 

26.9 

26.2 

2. East- West wind 

25.8 

25.5 

3. Density 

18.4 

22.8 

4. North-South wind 

5. Radial separation velocity 

20.8 

20.9 

(+0.02 m/s) 

9.2 

9.0 

6. Vertical wind 

1.5 

3.4 

7. C^: hypersonic-continuum (+4%) 3.1 

3.2 

8. C^: supersonic-continuum (±10%) .... 1.7 

1.9 

9. Initial mass (+0.5 kg) 

1.0 

1.0 


All other contributors < 0.5 km 
*3-o variance shown in parentheses 

persion analysis. Over 3200 random trajectories are sim- 
ulated to assure proper Gaussian or uniform distribu- 
tions for the 41 mission uncertainties identified. 

The statistical results from the 3200 Monte Carlo 
simulations are displayed in Figs. 8-15. Figures 8-10 
show the distribution of the total angle-of-attack at three 
discrete locations during the early phase of the mission: 
at atmospheric interface, in the transitional regime, and 
at peak heating. At atmospheric interface, the statistical 
mean total angle-of-attack of the 3200 Monte Carlo cas- 
es is 2.5 deg. The maximum otj observed is around 8 
deg (which is below the mission constraint of 10 deg). 
In the transitional regime, the total angle-of-attack does 
increase from atmospheric interface due to the free mo- 



Fig. 8. Distribution of total angle-of-attack at 
atmospheric interface resulting from over 3200 Monte 
Carlo simulation cases. 



Fig. 9. Distribution of total angle-of-attack in 
transitional regime resulting from over 3200 Monte 
Carlo simulation cases. 



ay, deg 


Fig. 10. Distribution of total angle-of-attack at peak 
heating resulting from over 3200 Monte Carlo 
simulation cases. 

lecular instability. The mean is 8. 1 deg, and the max- 
imum observed is 30.4 deg. The higher spin rate pre- 
vents continued growth in the total angle-of-attack; so 
that, by peak heating (where the SRC is stable), the mean 
damps to 2.5 deg as seen in Fig. 10. The maximum 
ctp observed at peak heating is 8.6 deg, which is below 
the mission constraint of 10°. 

Figures 11-13 show the distribution of the drogue 
and main parachute deployment conditions. The mean 
Mach number at drogue chute deployment is 1 .4, as seen 
in Fig. 1 1 . The minimum deployment Mach number en- 
countered is 1 .27, which is high enough to avoid the sig- 
nificant effects of the transonic dynamic instability. The 
corresponding mean total angle-of-attack at drogue chute 
deployment (see Fig. 12) is 3.6 deg. with a maximum 


7 




Mach number 


Fig . 11. Distribution of Mach number at drogue chute 
deployment resulting from over 3200 Monte Carlo 
simulation cases. 



Fig. 12. Distribution of total angle-of-attack at drogue 
chute deployment resulting from over 3200 Monte 
Carlo simulation cases . 

of 8.4 deg (well below the mission constraint of 30°). 
Figure 13 shows the distribution of the main parachute 
deployment altitude. The mean deployment altitude is 
3.1 km, with a minimum occurring at 2.21 km. 

Figures 14 and 15 show the resulting distributions 
in downrange and crossrange at landing for the 3200 
Monte Carlo cases, respectively. The minimum down- 
range is -49.4 km (short) from the nominal landing point, 
whereas the maximum downrange is 54.1 km (long). The 
maximum crossrange obtained is 18.9 km from the nom- 
inal landing point. The resulting 3 -a ellipse has a major 
axis of 83.5 km (-40.1 short, 42.2 long) in downrange 
and a minor axis of 29.2 km in crossrange. This foot- 
print is within the UTTR site; however, it is approach- 
ing the upper boundary limit. Within the assumptions of 



Fig. 13. Distribution of altitude at main chute deploy- 
ment resulting from over 3200 Monte Carlo simulation 
cases. 



Fig . 14. Downrange distribution at landing resulting 
from over 3200 Monte Carlo simulation cases . 



Fig. 15. Crossrange distribution at landing resulting 
from over 3200 Monte Carlo simulation cases. 
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namic instability is present which could induce a tum- 
bling motion prior to parachute deployment. 
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Fig. 16. Landing range dispersion resulting from over 
3200 Monte Carlo simulation cases. 


The solution selected to address the high altitude 
instability is to increase the SRC entry spin rate to 16 
rpm. The higher spin rate, although not eliminating the 
instability, increases the gyroscopic stability of the SRC 
thereby retarding the effects of the free molecular static 
instability. To address the transonic/subsonic instabili- 
ty, a drogue chute (having a diameter of 0. 828 m) is added 
and a deployment algorithm based on an accelerometer 
activated timer resulting in a Mach 1.4 deployment is 
defined. The drogue chute serves to stabilize the SRC 
until main parachute deployment, and is shown not to 
introduce an unacceptably large increase in the landing 
footprint. 


Table 5. Summary of Monte-Carlo Analysis 


Attitude dispersion 

Atmospheric interface 

Mean 

Min. 

Max. 

3-c 

a T , deg 

...2.5 

0.3 

8.0 

3.3 

Transitional regime ay, deg 

...8.1 

1.8 

30.4 

10.6 

Peak heating ay, deg 

Drogue chute deployment 

...2.5 

0.5 

8.6 

2.5 

ay, deg 

...3.6 

1.5 

8.4 

2.3 

Landing dispersion 

Landing downrange, km 

...0.4 

-49.4 

54.1 

42.4 (long) 
-40.1 (short) 

Landing crossrange, km 

..-0.2 

-18.7 

18.9 

14.6 

Total range, km 

... 1.2 

0.3 

54.5 

23.6 


the present analysis, a 99.7 percent probability exists that 
the SRC will land within this 3 -a footprint ellipse. Fig- 
ure 16 shows the landing location of all 3200 Monte 
Carlo cases. Table 5 summaries these results. 

SUMMARY 

A six-DOF analysis of the nominal Stardust SRC 
entry reveals that two aerodynamic instabilities result in 
unacceptable capsule dynamics during the descent. The 
first instability resides in the high altitude free molecu- 
lar regime, while the second appears toward the end of 
the entry in the transonic/subsonic flow regime. These 
instabilities, if not eliminated or at least suppressed, could 
lead to mission failure. In the free molecular regime, a 
static instability exists which produces large excursions 
in the angle-of-attack (approaching 70°) early in the 
nominal entry profile. If off-nominal attitudes or atti- 
tude rates exist at atmospheric interface, a backward entry 
is very possible. In the transonic/subsonic regime, a dy- 


For this mission, 41 potential uncertainties were 
identified which could impact the entry. Initial state vec- 
tor and atmospheric property (density, and North-South 
and East-West winds) uncertainties were found to pro- 
duce the greatest downrange dispersions on the order of 
20-25 km each. Uncertainties from bus separation and 
aerodynamics produced dispersion between 5-10 km 
each. All other uncertainties resulted in dispersion less 
than 1 km. 

A Monte Carlo analysis of over 3200 off-nominal 
trajectories shows that the SRC attitude near peak heat- 
ing and drogue chute deployment to be within Stardust 
program limits. The resulting 3 -G landing footprint ob- 
tained was 83.5 km (-40.1 short, 42.2 long) in down- 
range and 29.2 km in crossrange (which is within the 
Utah Test and Training Range boundaries). Within the 
assumptions of the present study, a 99.7 percent proba- 
bility exists that the Stardust SRC will land within this 
3-0 ellipse. 

The instabilities in the Stardust SRC were revealed 
too late in the design process. If identified earlier, these 
types of instabilities could be eliminated by considering 
alternative capsule configurations that avoid the need 
for corrective measures later in a program. Therefore, a 
case is made for including six-DOF entry trajectory anal- 
yses early in the conceptual design phase. 

Finally, the resolution of the Stardust SRC instabil- 
ities relies heavily on the validity of the Monte Carlo 
analysis more than any previous mission. Increased de- 
pendence on entry simulations for mission success plac- 
es considerable importance on selecting appropriate un- 
certainties. As confidence increases in the analysis ac- 
curacy, cheaper and/or higher performance entry sys- 
tems can be selected for future missions. 
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